Three Dimensional Magnetic Correlations in Multiferroic LuFe2 04 
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We present single crystal neutron diffraction measurements on multiferroic LuFe2 04. Magnetic 
reflections are observed below transitions at 240 and 175 K indicating that the magnetic interactions 
in LuFe204 are 3-dimensional (3D) in character. The magnetic structure is refined as a ferrimagnetic 
spin configuration below the 240 K transition. Below 175 K a significant broadening of the magnetic 
peaks is observed along with the build up of a diffuse component to the magnetic scattering. 

PACS numbers: : 77.84.-s, 75.30.Kz, 28.20.Cz, 25.40.Dn 
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Materials that offer the possibility of simultaneously 
controlling magnetic and electric degrees of freedom are 
the subject of intense interest [H . Recently, multiferroic 
materials have been identified that show large coupling 
between electric and magnetic degrees of freedom. Ferro- 
electricity driven by either magnetic or charge ordering 
appears to be the origin of the large coupling and, hence, 
understanding the underlying electronic interactions is 
crucial for further insight into multiferroicity 

LuFe204 has attracted attention as a novel ferroelec- 
tric material where ferroelectricity is driven by the elec- 
tronic process of charge ordering of Fe^+ and Fe'^"'' ions 
and for indications of coupling between electronic and 
magnetic degrees of freedom 0, S S S Si- LuFe204 is 
a member of the i?Fe204 {R — rare earth element) fam- 
ily, the physical properties of which depend strongly on 
oxygen stoichiometry. For example, nearly stoichiometric 
YFe204 exhibits three-dimensional (3D) magnetic order 
while oxygen deficient YFe2 04 exhibits two-dimensional 
(2D) magnetic order Q. LuFe204 exhibits multiple 
phase transitions. 2D charge correlations are observed 
below 500 K, while below 320 K 3D charge order is es- 
tablished, roughly coinciding with the onset of ferroelec- 
tricity [3, 01 . Magnetic order appears below 240 K and 
2D ferrimagnetic order has been suggested by neutron 
scattering studies Q. However, strong sample dependent 
behavior observed in other members of i?Fe204 [7| sug- 
gests that unraveling the interesting behavior of LuFe2 04 
requires paying due attention to sample quality. 

In this letter we present extensive neutron diffraction 
measurements from 20 to 300 K on high quality single 
crystals of LuFe204. We report several new findings that 
provide information about the underlying magnetic in- 
teractions. First, our measurements indicate that below 
240 K 3D magnetic correlations exist with magnetic in- 
tensity appearing at (1/3 1/3 L) where L may take on 
integer and 1/2-integer values. The magnetic structure 
is refined with a ferrimagnetic spin configuration with 
a propagation vector of (1/3 1/3 0). The magnetic in- 
tensity appearing on peaks where L is a 1/2-integer is a 



consequence of the charge ordering at ~320 K. In addi- 
tion, evidence is presented for a second transition at 175 
K with significant changes in magnetic peak intensities 
and broadening of many reflections. 

Single crystals of LuFe204 were grown by floating- 
zone-melting, using an oxygen partial pressure tune d by 
a CO/CO2 mixture to control oxygen stoichiometry [lO|. 
For CO/CO2 ratio close to 2.7 the magnetic behavior 
exhibits two sharp magnetic transitions in contrast to 
previous single crystal magnetization measurements Q 
which show only a single transition. The magnetic be- 
havior is qualitatively similar to that of stoichiometric 
YFe204 where previous work with slightly varying oxy- 
gen concentration showed that stoichiometric samples are 
characterized by much sharper magnetic transitions and 
reduced residual low temperature (T) susceptibility 11 1 . 
We take the sharpness of the observed transitions and 
reduced residual low temperature susceptibility, in light 
of previous measurements on YFe204, as strong evidence 
that these LuFe204 crystals have nearly ideal oxygen sto- 
ichiometry and are extremely homogeneous |12l |. Two 
crystals of the growth batch in which the sharpest tran- 
sitions were observed, denoted SI and S2, were selected 
for neutron diffraction measurements. S2 is slightly more 
homogeneous as judged from magnetization curves. 

Neutron diffraction measurements were performed on 
the N5 triple-axis spectrometer (TAS) at the Canadian 
Neutron Beam Center at Chalk River Laboratories and 
the HBl TAS at the High Flux Isotope Reactor at Oak 
Ridge National Laboratory. Neutrons with incident ener- 
gies of 14.56 meV (N5) and 30.5 meV (HBl) were selected 
and horizontal coUimations of 30'-36'-16.4'-66' (N5) and 
48'-40'-40'-70' (HBl) were used. Pyrolytic graphite (PC) 
(002) was used for monochromator and analyzer and PG 
filters were placed in the scattered beam to suppress 
higher order contamination. 

Figure 1 shows the integrated intensity vs. T for the 
(1/3 1/3 0) magnetic peak providing evidence for two 
phase transitions, one at 240 K (T/v) and another at 175 
K [Tl) [3]. The presence of two transitions is corrob- 
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FIG. 1: (Color online) Integrated Intensity for the (1/3 1/3 0) 
magnetic peak measured in the (HKO) plane. The error bars 
are smaller than the point size. The inset shows field-cooled, 
c-axis magnetization data at 0.03 (red), 0.07 (black), and 0.1 
Tesla (blue). Arrows indicate whether the measurement was 
made upon warming or cooling. 
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FIG. 2: (Color online) (a) and (c) show scans along (1/3 1/3 
L) for SI and S2. The solid blue lines are described in the 
text. The error bars in (c) have been omitted for clarity but 
are comparable to those displayed in (a). The arrows indicate 
peaks contaminated by Aluminum background scattering, (b) 
Displays data showing the appearance below Tl of a new set 
of satellites indexed by (l/3±(5 l/3±5 3L/2) with 5 ~0.027. 



orated by the c-axis magnetization (inset Fig. 1) of a 
crystal from, the same batch as the neutron scattering 
samples. The sharpness of in both neutron and mag- 
netization data suggests that Tl is of 1st order, which is 
confirmed by hysteresis in the magnetization. 

Figure 2 (a) and (c) displays scans along (1/3 1/3 L) at 
several temperatures. At 280 K peaks at large values of 
L are readily visible with 1/2-integer indices (Fig. 2(c)). 



As pointed out previously 0], these peaks are consistent 
with a -\/3 X -\/3x 2 unit cell containing 36 Fe atoms. 
Measurements at larger values of momentum transfer, Q, 
indicate that the peaks are not magnetic and, in accord 
with previous work, they are attributed to the onset of 
3D charge order at 320 K Q. 

Cooling below 240 K, new intensity appears at integer 
values of L. This intensity is strongest at small L and for 
L > 16 it has diminished to the point where it is difficult 
to discern from the background. Measurements at (2/3 
2/3 L) and (4/3 4/3 L) reinforce the conclusion that the 
intensity on integer L positions only occurs for sufficiently 
small values of Q. Similarly, the intensity at 1/2-integer 
positions become enhanced at small Q, but not at larger 
Q where the intensity is essentially unchanged from that 
observed at 280 K. Such Q-dependence is expected for 
scattering from the magnetic moments of Fe^+ and Fe'^"'" 
indicating that peaks indexed by integer values of L orig- 
inate principally from an ordered magnetic sublattice. 
This also indicates that the changes of the L 1/2-integer 
peaks are predominantly due to magnetic order rather 
than modifications of the charge ordering configuration. 
Thus, the neutron scattering data demonstrates 3D mag- 
netic correlations below 240 K. However, we do note that 
the magnetic reflections are not resolution limited along 
L suggesting a finite correlation length. 

To allow for quantitative comparison with models for 
the spin configuration, a large number of reflections were 
measured at 220 K by scanning along the (1/3 1/3 L), 
(2/3 2/3 L), and (4/3 4/3 L) directions. To solve for the 
magnetic structure, representational analysis was per- 
formed to consider those magnetic structures which are 
syinmetry allowed from the parent i?3m space group 
[3, [iBl ■ This analysis assumes that the onset of charge 
order at 320 K does not signiflcantly affect the symmetry- 
allowed magnetic order. Initially, a magnetic propaga- 
tion vector of (1/3 1/3 1/2) was considered and symme- 
try analysis yielded two possible irreducible representa- 
tions. However, both possibilities yielded L 1/2-integer 
reflections several orders of magnitude more intense than 
the L-integer reflections. This is in contrast to observa- 
tion where the 1/2-integer peaks are comparable to or 
weaker than the integer reflections (see Fig. 2 (c)). Con- 
sequently, we concluded that the magnetic structure is 
described by the ordering wavevector (1/3 1/3 0), the 
presence of 1/2-integer reflections occurring as a result 
of the charge ordering which decorates the lattice with 
differing magnetic moment on Fe^+ and Fe^+ sites with 
a periodicity of (1/3 1/3 1/2). Representational analysis 
with the (1/31/30) wavevector again yielded two allowed 
irreducible representations. For spins pointing along the 
c-axis, as suggested by the magnetization measurements, 
these representations correspond to ferromagnetic (FM) 
or antiferromagnetic (AFM) alignment of the two spins 
of the primitive basis (see Fig. 3). The AFM case can 
be ruled out immediately as the magnetic structure, in- 
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FIG. 3: (Color online) Magnetic structure of LuFe204 at 
r=220 K. The three magnetic domains correspond to the 
three symmetry equivalent magnetic propagation vectors (1/3 
1/3 0), (1/3 -2/3 0), and (-2/3 1/3 0) labeled as Dl, D2, and 
D3 respectively. Bl and B2 denote the primitive basis ((0, 0, 
0.22) and (0, 0, 0.78)). 
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FIG. 4: (Color online) Peak intensity for scans along (1/3 1/3 
L), (2/3 2/3 L) and (4/3 4/3 L) corrected as described in the 
text. The solid red line represents the result of the model 
calculation described in the text. 

eluding symmetry equivalent wavevectors, does not yield 
intensity at the (1/3 1/3 0) position. The FM coupling 
between spins in the basis, the only remaining symmetry 
allowed possibility, results in a ferrimagnetic structure as 
shown in Fig. 3 for the symmetry equivalent propagation 
vectors (1/3 1/3 0), (-2/3 1/3 0), and (1/3 -2/3 0). 
Determining the agreement between the model and ex- 



periment was complicated by reflections that were not 
resolution limited along the c-axis. To account for this, 
the TAS resolution function was fully simulated [l3| and 
the data from several reflections were fit assuming a 
Lorentzian lineshape along the c-axis with resolution lim- 
ited peaks in the hexagonal basal plane. This nicely re- 
produced the lineshape and yielded a Lorentzian HWHM 
(A) of 0.0257 r.l.u. corresponding to a correlation length 
( ^ = 1/A) of 156 A. The resolution correction was then 
obtained by convolving this lineshape with the resolu- 
tion function for all measured reflections. For simplicity, 
the magnetic moments were fixed to be the same on all 
sites. This assumption is justified given the charge or- 
der which occurs at 320 K and the domain population of 
this charge ordered state. To accurately reproduce the 
1/2-integer L peaks, the correct charge ordering scheme 
would need to be included, but it is not necessary to de- 
scribe the integer-L reflections. Therefore we proceeded 
to fit 58 observed integer-L reflections with a model us- 
ing only 4 parameters, two domain population factors, 
an overall scale factor, and a Debye- Waller factor. The 
model agrees very well with the data (reduced =1.39). 

To extract meaningful domain populations, it was rec- 
ognized that the Lorentzian width varied slightly between 
the three domains yielding correlation lengths of 155 A, 
130 A, and 160 A for domains Dl, D2, and D3 (see Fig. 
3). The resolution correction was modified to account 
for this and the resulting data corrected for resolution, 
magnetic form factor, and spin polarization factor are 
plotted in Fig. 4 along with the calculated intensities. 
The "^as improved slightly to 1.37 and the resulting 
domain population ratios are 0.85:0.57:1. Comparable 
domain population ratios are observed for SI, a sample 
taken from the same growth. Whether or not the domain 
populations are an intrinsic property of LuFe2 04 or are 
sample dependent is, at present, unclear. Domain Dl 
yields a symmetric pattern in L with peaks described by 
L=3n while D2 (D3) produce peaks at L=3n-|-1 (3n-l) for 
H=l/3 and 4/3, the order reversed for H=2/3. The ratio 
of D2:D3 (0.57:1) is responsible for the asymmetry in the 
diffraction pattern for ±L. The proposed magnetic struc- 
ture shown in Figure 3 is ferrimagnetic with an excess of 
1/3 of the spins pointing along the c-axis. This results 
in a saturation ordered moment of 1/3 x 4.5 /is/Fe x 
2 Fe/FU or 3 /xs/FU. This value is very close to the 
saturation magnetization in higher fields {H > 2 T) in 
our crystals (not shown) as well as to the value reported 
in Ref. Q further corroborating the magnetic structure. 
The smaller magnetic moment observed in the magneti- 
zation at if < 2 T is a consequence of the formation of 
domains and will be discussed in more detail elsewhere. 

Below Ti an additional component to the scattering 
builds up which is extremely broad along (1/3 1/3 L) but 
sharp along (HHO) (Fig 2(a) and (c)). This diffuse scat- 
tering appears to be magnetic in origin as emphasized 
by the solid blue lines in Fig. 2 (a) and (c) which are 
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proportional to the product of the form factor squared, a 
Debye-WaUer factor, and the polarization factor (l+Q^). 
The diffuse scattering is stronger in SI which, as judged 
by magnetization measurements, is slightly less homoge- 
nous than S2 suggesting that the presence of 2D magnetic 
short range order is a property of samples that are not 
sufhciently homogeneous. This bears some similarity to 
YFe2 04 where the dimensionality of the magnetic inter- 
actions (2D or 3D) depends on the oxygen stoichiometry 
• Figure 2 also shows that below Tl profound changes 
occur in the magnetic peaks along (1/3 1/3 L). The in- 
tensity along (1/3 1/3 L) for magnetic reflections changes 
rather dramatically with some peaks becoming more in- 
tense (e.g. (1/3 1/3 1)) and some peaks becoming less 
intense (e.g. (1/31/3 0)). Thus 3D magnetic correlations 
persist below albeit with a shorter correlation length 
than found for T/v. Finally, we note that scans along 
(110) have revealed the existence of a new set of satellite 
peaks of unknown origin indexed as (l/3±(5 l/3±(5 3L/2) 
where S ~0.027 (see figure 2 (b)) below Tl. 

The increase in linewidth on many peaks below 
is consistent with the introduction of stacking faults at 
a structural phase transition |l3|. Evidence that Tl 
involves a structural component is provided by an ex- 
tinction related increase in intensity of strong structural 
Bragg peaks as well as broadening of those peaks be- 
low 175 K. CuKa powder x-ray diffraction as well as the 
single crystal neutron data shows a contraction (expan- 
sion) of the a (c) lattice constant from 350 K to 100 K 
but does not reveal any sign of a structural distortion. 
Stoichiometric YFe2 04 provides an example of a related 
system which exhibits phase transitions with a structural 
component at similar temperatures to the two transitions 
found in LuFe204 below 300 K ^. 

Although the linewidth of the peaks below Tl and the 
broad diffuse component to the scattering makes a full 
solution of the low T magnetic structure very difficult, 
some general conclusions are possible. The 3D magnetic 
correlations in LuFe2 04 are intrinsically sensitive to dis- 
ruptions along the c-axis as the superexchangc path be- 
tween Fe-0 bilayers must pass through not one but two 
oxygen ions and thus the driving force for the changes in 
the magnetic structure below Tl may be related to the 
introduction of stacking faults as discussed above. The 
stacking arrangement would result in local magnetic or- 
der which deviates from the ferrimagnetic state and may 
mix in states which are locally AFM. If we consider the 
symmetry allowed AFM state, the most intense reflection 
is the (1/3 1/3 ±1) reflection, the most strongly enhanced 
and broadened reflection seen experimentally in the low 
T state. Additionally, the (1/3 1/3 0) peak is absent in 
the AFM structure, qualitatively consistent with the ob- 
served strong reduction in the intensity of this reflection 
seen at low temperatures. Furthermore, local AFM order 
would result in decreased c-axis magnetization, leading 



naturally to the similar temperature dependence of the 
magnetization and the (1/3 1/3 0) intensity (Fig. 1). A 
detailed quantitative analysis of the neutron scattering 
and magnetization data below 175 K, beyond the scope 
of the present work, should shed further light on the com- 
plex magnetic behavior in LuFe204. 

In conclusion, we show that that LuFe204 has two 
transitions below 300 K. Both transitions involve a 3D 
magnetically correlated structure with a finite correla- 
tion length along the c-axis. Whether the correlation 
length is an intrinsic property or is the result of disor- 
der, most likely oxygen stoichiometry, is yet to be elu- 
cidated. Below T/v a ferrimagnetic spin configuration is 
found with a magnetic propagation vector of (1/3 1/3 0) 
with magnetic intensity occurring at (1/3 1/3 L) where L 
is 1/2-integer arising due to the charge ordering at 320 K. 
Theoretical models taking into account the 3D nature of 
the magnetic interactions as well as the sequence of mag- 
netic phase transitions described above should provide 
insight into the multiferroic behavior of LuFe204. 
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